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Notati on 


Roman  Letter  Symbol s 

a  Radius  of  dipole 

S’  Magnetic  vector  potential 

D  Distance  from  cable  to  intruder 

f  Electric  field 

f  Frequency 

F  (  )  Pulse  expansion  functions 

G(  )  Free  space  scalar  Green's  function 
h0  Transverse  wave  number 

H  Scalar  Magnetic  field,  e  component 

0 

)  Hankel  function  of  the  first  kind 

T  Line  current 

7  Surface  current  density 

k0  Free  space  propagation  constant 

Kj(  )  Modified  Bessel  function  of  the  second  kind,  order 
one 

L  Length  of  dipole 

N  Number  of  sub-dipoles 

t  Time 

v  Phase  velocity 

X  Distance  from  cable  to  receiving  antenna  (12 
wavelengths) 

Azn 


o 


Length  of  each  sub-dipole 


Notation 
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Greek  Letter  Symbo I s 
a  Angle  of  incidence  on  dipole 

8  Propagation  constant 

Yn  Angle  of  incidence  on  receive  antenna,  from  dipol 

e  Relative  permittivity 

e„  Free  space  permittivity 

n0  Free  space  characteristic  impedence 

6  Angle  of  departure  from  dipole 

X  Wavelength 

p  Relative  permeability 

y0  Free  space  permeability 

p  Distance 

o  Conductivity 

<j>  Angle  of  incidence  on  receive  antenna,  from  cable 

<j>  Scalar  potent  ' 1 

oj  Angular  frequency 
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Abstract 

The  free  space  interaction  of  the  near  zone 
electromagnetic  fields  external  to  a  leaky  coaxial  cable  and 
a  moving  dipole  is  investigated.  The  method  of  moments  and 
the  equivalence  principle  are  used  to  find  the  currents  on 
the  dipole  and  the  electromagnetic  fields  scattered  from  the 
dipole.  Finally,  the  overall  interaction  between  these  two 
fields  is  determined  by  superposition.  The  theoretically 
derived  results  are  then  compared  with  experimental 
measurements . 

The  findings  of  this  report  indicate  that  exact  results 
are  not  obtained  with  a  free  space  analysis;  therefore, 
further  work  involving  electronic  intruder  detection 
analysis  should  include  the  effects  of  the  earth-air 
interface . 


vi  i 


FREE  SPACE  ANALYSIS  OF  THE  SCATTERED 


ELECTROMAGNETIC  FIELDS  AS  A  DIPOLE 
APPROACHES  A  LEAKY  COAXIAL  CABLE 


I .  INTRODUCTION 


1 • 1  Background 

The  USAF,  DOD  and  many  other  government  agencies  are 
investigating  the  use  of  electronic  intruder  detection 
systems  for  protection  of  critical  permanent  and  semi¬ 
permanent  areas,  such  as  MX  sites,  nuclear  power  plants, 
forward  located  air  strips,  and  other  high  value  resources. 
One  such  system  is  the  Single  Wire  Individual  Resource 
Protection  Sensor  (SWIRPS).  The  SWIRPS  system  consists  of  a 
single  leaky  coaxial  cable,  which  encircles  the  resource  to 
be  protected,  and  a  central ly- 1 ocated  monopole  receiving 
antenna.  The  leaky  coaxial  cable  allows  a  small  portion  of 
the  RF  energy  travelling  inside  the  cable  to  "leak"  out, 
either  through  gaps  in  a  loose  braid  or  through  periodic 
slots  cut  in  a  solid  outer  conductor,  and  propagate  along 
the  outside  of  the  cable.  As  an  intruder  approaches  the 
leaky  coaxial  cable,  which  acts  as  a  waveguiding  structure, 
the  Electromagnetic  (EM)  fields  are  disturbed,  thereby 
changing  the  received  signal  at  the  monopole.  When  the 


change  in  the  received  signal  is  sufficiently  large  a 
detection  is  declared;  therefore,  it  is  the  interaction  of 
an  intruder  with  the  external  EM  field  that  allows 
detection. 

1 . 2  Problem 

The  intent  of  this  thesis  is  to  study  the  interaction 
between  a  radially  approaching  intruder  and  the  EM  fields 
created  by  the  leaky  coaxial  cable,  so  that  the  resulting  EM 
fields  at  the  receiving  antenna  can  be  determined.  In  order 
to  scope  down  the  problem  to  a  workable  and  more  general 
case,  only  the  plane  containing  the  intruder,  the  coaxial 
cable,  and  the  receiving  antenna  will  be  examined. 

Before  any  analysis  of  the  interaction  between  the  EM 
fields  and  the  intruder  can  be  made,  the  propagation 
characteristics  of  a  leaky  coaxial  cable  must  be  known.  A 
solution  for  lossless  leaky  coaxial  cables  in  free  space  has 
been  obtained  by  Fernandes  (Ref  4).  His  solution  is  based 
on  the  concept  of  surface  transfer  impedence  and  shows  the 
existance  of  two  propagating  modes,  the  coaxial  and  the 
'Goubau'  modes.  Fernandes  indicates  that  the  coaxial  mode 
exists  when  the  cable  is  excited  in  the  normal  coaxial 
manner;  if  however,  the  cable  is  excited  as  a  solid  wire, 
the  'Goubau'  mode  exists.  The  coaxial  mode  has  the  majority 
of  its  energy  propagating  inside  the  cable,  is  little 
affected  by  the  outside  medium,  and  thus  is  of  minimal  help 
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in  intruder  detection.  The  'Goubau'  (or  surface  attached) 
mode  however,  has  a  bound,  non- radi ati ng ,  radially 
evanescent  field  with  the  majority  of  its  energy  propagating 
outside  the  cable,  and  is  therefore  very  dependent  upon  the 
outside  environment.  The  radial  decay,  or  attenuation,  of 
the  'Goubau'  mode,  though  highly  dependent  on  the  frequency 
of  operation  and  the  inside  dielectric  permittivity 
(Ref  4:259),  is  also  affected  by  the  outside  environmental 
conditions,  such  as  snow,  rain,  ice,  etc  (Ref  13:135),  which 
are  changing  conditions  that  cannot  be  designed  in  or 
permanently  fixed.  Future  work  should  investigate  the 
effect  of  the  cable  located  on  or  in  an  earth  model,  which 
could  also  be  used  to  predict  any  detection  degradation  due 
to  changing  environmental  conditions. 

Some  initial  analysis  of  the  i ntruder/1 eaky-cabl e 
interaction  have  been  performed  (Ref  8:7-17,  Ref  10); 
wherein  the  intruder  was  simply  modeled  as  a  scattering 
point  source,  therefore,  no  effects  due  to  size  or 
orientation  could  be  seen.  As  an  obvious  continuation  of 
the  previous  works,  this  thesis  will  model  the  intruder  as  a 
dipole  thereby  making  it  possible  to  determine  the  effects 
of  intruder  orientation.  The  vertical  or  standing  intruder 
is  the  only  orientation  considered  herein.  The  effects  of  a 
non-plane  wave  incident  upon  a  dipole  must  be  developed 
since  the  intruder-dipole  model  is  located  in  the  near  field 
region  of  the  leaky  coaxial  cable. 
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The  SWIRPS  system  being  considered  has  the  cable 
located  on  the  surface  of  the  earth;  thus,  the  leaky  coaxial 
cable,  the  receiving  antenna,  and  the  intruder  are  on  the 
same  base  plane.  Also,  since  experimental  measurements 
indicate  that  the  zone  of  detection  is  less  than  a 
wavelength  (Ref  11:18),  the  intruder  will  only  be  considered 
to  move  within  + two  wavelengths  of  the  coaxial  cable. 

1 • 3  Assumptions 

An  analysis  of  this  effort  will  be  made  making  the 
following  assumptions: 

1.  All  analysis  will  be  performed  in  free  space.  This 

assumption  is  being  made  since  the  only  form  available  for 

the  external  EM  fields  from  a  leaky  coaxial  cable  are  for 

the  cable  located  in  free  space,  or  for  the  cable  located 
within  a  mine  shaft.  It  was  decided  that  the  free  space 
analysis  would  be  more  appropriate  to  this  study.  In 
addition,  the  skin  depth  for  dry  clay  soil  at  10  MHz  (Ref 
16:314)  was  calculated  to  be  in  the  range  of  85  meters,  thus 
indicating  minimal  effect  of  the  earth-air  interface  since 
the  radial  electric  field  normally  decays  at  a  much  faster 
rate . 

2.  The  leaky  coaxial  cable  is  considered  to  be 

lossless,  infinitely  long,  and  straight. 

3.  The  intruder  will  be  modeled  as  a  perfectly 

conducting  thin  cylinder  (approaching  a  dipole).  The 
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intruder  dipole  model  will  be  0.355  wavelengths  in  length 
with  the  radius  being  much  less  than  a  wavelength;  since  the 
operating  frequency  is  60MHz  (A=5  meters),  this  corresponds 
to  a  5'-10"  intruder. 

1.4  Approach  and  Presentation . 

The  theory  is  developed  in  Chapter  II  using  the 
following  approach  to  determine  the  interaction  of  a  moving 
intruder  with  the  EM  fields  created  by  the  leaky  coaxial 
cable.  First,  the  field  induced  currents  on  the  intruder 
must  be  determined;  this  will  be  accomplished  via  the  method 
of  moments.  The  method  of  moments  being  a  technique  for 
solving  an  integral  equation,  which  can  be  readily 
implemented  on  a  computer.  Once  an  algorithm  is  developed 
to  determine  the  field  induced  currents  on  the  intruder,  the 
electric  field  at  the  receiving  antenna  can  be  determined  by 
summing  the  source  field  and  the  intruder  scattered  field. 
Finally,  the  interaction  of  the  leaky  coaxial  EM  fields  and 
the  intruder  (with  position)  can  then  be  obtained.  Chapter 
II  also  contains  sample  numerical  results  giving  the 
incident  field  along  the  dipole  and  the  resulting  currents 
induced  on  the  dipole  for  two  selected  positions. 

Chapter  III  presents  the  results  obtained  from  the 
developed  computer  program  (Appendix  A)  using  the  approach 
given  in  Chapter  II.  The  received  power  at  the  receiving 
antenna  as  a  function  of  intruder  distance  from  the  cable  is 
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given  for  two  ambient  power  levels.  These  results  are  then 
compared  to  those  calculated  and  experimentally  obtained  by 
Poi rier  (Ref  8) . 

The  conclusions  and  recommendations  are  presented  in 
Chapter  IV. 


II.  DEVELOPMENT  OF  THEORY 


Loose  braid  or  'leaky1  coaxial  cables  have  been  studied 
primarily  for  the  purposes  of  communications  within  mine 
shafts.  The  theory  developed  in  these  cases  is  not 

applicable  to  this  study  because  the  mine  shaft  problem 
imposes  a  different  boundary  condition  on  the  field 

development;  that  is,  for  the  case  of  the  mine  shaft,  the 

fields  are  such  that  they  must  go  to  zero  at  the  tunnel 

walls,  as  opposed  to  going  to  zero  at  infinity. 

The  propagation  characteristics  for  the  special  case  of 
a  loose  braid  coaxial  cable  in  free  space  has  been 
determined  by  Fernandes.  His  solution,  which  is  based  on 
the  surface  transfer  impedence  concept,  shows  the  existence 
of  two  propagating  modes,  the  coaxial  mode,  and  the  'Goubau' 
mode.  Fernandes  describes  the  fields  external  to  the  leaky 
cable,  for  either  mode,  as 

E_,  =  A  o  H i 1 jh0p)exp(-jgz)exp( jut)  (2.1) 

E  =  Aor —  H  j 1  ^  ( j h o p ) e xp ( - j Bz ) e xp ( juit )  (2.2) 

P  no 

k  o 2  ( j  \ 

H0  =  A ° lo~u~qTi~  H*  (ohop)exp(-jez)exp( jwt) 


(2.3) 


where  )  are  hankel  functions  of  the  first  kind  for 

order  zero  and  one,  0  is  the  propagation  constant  ^which 
is  determined  by  the  mode  selected),  k0  is  the  free  space 
propagation  constant,  h0  is  the  transverse  wave  number 
defined  as  ho2=02-ko2,  and  A0  is  a  constant  assumed  to  be 
unity  hereafter  iRef  4:256). 

Since  the  problem  analyzed  in  this  study  is  only  for  a 
particular  cut  of  an  infinitely  long,  lossless,  and  straight 
cable,  the  phase  variation  arising  from  the  e”J^z  term  is 
constant  for  any  particular  'z‘,  and  thus  can  be  ignored. 
Also,  realize  that  for  the  case  of  a  lossless  line,  as 
considered  herein,  h0  is  always  a  real  number,  and 
therefore  by  using  the  relations  between  solutions  for 
Hankel  functions  and  the  Modified  Bessel  functions  of  tne 
second  kind,  (Ref  1:375),  the  expressions  for  Ep  and  HQ 
become 

ep  ■  F7  j-  I  <2-4> 


He 


k02 

u>Uoh  o 


2 

IT 


K i (h0p) 


(2.5) 


Experimental  measurements  indicate  the  existence  of  two 
interfering  waves,  propagating  in  the  z  direction,  external 


to  the  leaky  cable  (Ref  9);  this  implies  that  a  form  of  the 


'Goubau'  mode  is  excited  by  the  coaxial  mode  fields  that 
have  'leaked'  to  the  outside.  Based  on  the  above 
assumption,  the  field  structure  used  in  this  study  is  of  the 
form 


Ep  =  E p ( 3 x )  +  C  Ep( 32 )  (2.6) 

where  6i  is  the  propagation  constant  for  the  coaxial  mode, 
g2  is  the  propagation  constant  for  the  'Goubau'  mode,  and 
C  is  a  constant  assumed  to  be  unity  hereafter.  The  values 
for  gj  and  g2  are  calculated  using  the  relationship 
,  where  v  is  the  phase  velocity  and  f  is  the 
frequency  of  operation.  Experimental  measurements  of  Bi 
and  62  have  been  made  at  75  MHz  (Ref  8:23);  thus,  the 
values  of  Bi  and  62  at  60  MHz  can  be  calculated  as 
81  =1.72  rads/m  and  62=1.32  rads/m 

2 . 1  Non-Plane  Wave  Incidence  on  a  Dipole 

The  moving  dipole,  as  shown  in  Figure  2.1,  will 
interact  with  the  vertical  component  of  the  electric  field 
described  by  equation  2.6.  The  close  proximity  of  the 
dipole  to  the  leaky  coaxial  cable  results  in  a  significant 
variation  of  the  electric  field  across  the  length  of  the 
dipole.  In  order  to  determine  the  excitation  of  the  dipole, 
we  must  consider  the  dipole  as  being  broken  into  N  equal 
sub-dipoles,  or  segments,  where  the  length  of  each  segment. 


L/N,  is  such  that  the  electric  field  variation  across  it  is 
minimal.  Also,  assume  that  the  electric  field  present  at 
the  center  of  each  segment  is  also  the  electric  field  that 
exists  over  the  entire  segment.  The  dipole  is  now  being 
viewed  as  an  array  of  N  dipoles,  each  being  excited  by  an 
incident  electric  field  which  can  be  described  as 

eJ  =  [Ep<0i )  +  Ep(B2)]sin(an)  (2.7) 

where  E  j!,  is  the  vertical  component  of  the  incident 

j.  i_ 

electric  field  at  the  nL  element  and  a  is  the  angle 

n 

of  incidence  at  the  n*"*1  element,  n  being  1  at  the  bottom 
and  N  at  the  top  of  the  dipole.  The  angle  of 
incidence,  ap  ,  is  a  function  of  the  dipoles  horizontal 
distance  from  the  cable,  D,  as  well  as  the  segment  being 
excited,  n  ,  and  is  given  by 


a 


n 


arc  tan 


(2.8) 


Thus,  the  amount  of  coupling  at  each  sub-dipole  or  segment 
is  seen  to  be  dependent  upon  the  location  of  the  dipole  and 
also  the  location  of  the  segment  being  excited.  Figure  2.2 
is  a  plot  of  the  magnitude  of  the  incident  electric  field  by 
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Figure  2.2  Incident  Electric  Field  Variation  Across  the 
Dipole  for  Two  Intruder  Distances,  D  =  0.05  and  D=  1 .25 
Wavelengths.  (Note:  The  units  of  the  electric  field  are 

vol ts/wavel ength ) . 
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segment  for  two  distances  from  the  coaxial  cable,  D=0.05  and 
D=1.25  wavelengths.  For  small  distances,  as  shown  in  the 
D=0.05  curve,  the  variation  of  the  sine  function  overshadows 
the  Bessel  function  variation. 

Now,  in  order  to  determine  the  induced  dipole  current, 
the  method  of  moments  will  be  used;  it  is  therefore 
convenient  that  the  incident  electric  field  has  already  been 
determined  at  each  of  the  N  segments. 

2 . 2  Development  of  the  I nduced  Dipole  Current 

The  development  presented  in  this  section  is  a 
redevelopment  of  the  pertinent  sections  of  chapter  7,  Ref 
14.  The  induced  dipole  current  gives  rise  to  a  scattered 
electric  field,  Fs  ;  so,  in  order  to  obtain  an  expression 
for  Fs  ,  the  induced  dipole  current  will  be  determined 
using  the  method  of  moments.  The  method  of  moments  is  a 
computer  oriented  technique  which  enables  one  to  solve  an 
integral  equation,  in  this  case,  Pocklington's  Integral 
Equation,  in  terms  of  the  desired  dipole  current.  The  first 
step  is  to  approximate  the  integral  equation  by  a  set  of 
linear  algebraic  equations  with  the  source  induced  current 
on  the  dipole  being  the  unknown.  The  dipole  is  then  divided 
into  N  segments  with  each  segment  given  an  unknown  current 
function;  for  this  study,  pulse  function  currents  of  unknown 
amplitude  will  be  used  as  the  unknown  current  functions 
because  of  the  greatly  simplified  integration  associated 
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with  their  use.  The  problem  is  now  expressed  in  the  form 
V=ZI  ,  where  V  represents  the  incident  electric 
field,  E"1  ;  Z  is  a  form  of  an  impedence  matrix;  and,  I  is 
the  unknown  current  distribution.  An  expression  for  the 
unknown  current  distribution,  in  terms  of  the  source  or 
incident  field,  can  now  be  obtained  using  a  matrix  inversion 
technique.  Since  pulse  function  currents  are  used,  the 
unknown  current  distribution  is  only  approximated;  but  if 
the  width  of  the  pulse  functions  is  much  less  than  a 
wavelength,  a  good  approximation  of  the  dipole  current 
distribution  is  obtained  by  curve  fitting  the  values  of  the 
pulses.  Also,  with  this  form  of  excitation,  the  induced 
current  had  a  linear  6°  phase  variation  across  the  dipole. 
Figure  2.3  is  a  plot  of  the  induced  current  by  segment  for 
the  two  incident  electric  fields  plotted  in  Figure  2.2;  note 
that  while  the  induced  current  distribution  maintains  the 
same  basic  shape,  the  maximum  level  is  skewed  towards  the 
segment  which  had  maximum  incident  electric  field. 

2.2.1  The  Pocklington  Integral.  Assume  a  dipole  of 
length  L  and  radius  a  is  located  in  free  space,  as  shown 
in  Figure  2.4.  Also,  assume  that  the  radius  of  the  dipole 
is  much  less  than  the  free  space  wavelength  and  the  dipole 
length  L  ;  it  is  then  reasonable  to  assume  that  the 
induced  dipole  current  has  only  a  z  component.  In  this 
case,  the  Lorentz  Condition  v?f=  -jcoyoe0$  will  reduce  to 
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Figure  2.3  Induced  Current  Distribution  for  Two  Intruder 
Distances,  D=0.05  and  D=1.25  Wavelengths.  (Note:  The  units 
of  the  induced  current  are  amps/wavelength). 
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Figure  2.4  Dipole  in  Free  Space. 


aAz 

—7“  =  -jwpoEo^ 


where  <t>  is  the  scalar  potential  and  is  the 

directed  component  of  the  magnetic  vector  potential. 

Now,  writing  Maxwell's  Equation  in  terms  of 
magnetic  potential  A  as 


E  = 


-juft  + 


j  <>M*  11  C  0 


( 


and,  since  ft  =  \  a  ,  equation  2.10  can  be  reduced  using 
the  relationship  of  equation  2.9  to 


E 


Z 


-jo)Az 


34 

3z' 


From  equation  2.9  solve  for 


84 
8  z 


as 


(2.11) 


ii  -  .  _  .  1  a2Az  (2.12) 

8  z  jfcJpoeo  8ZZ 


Now,  substituting  equation  2.12  into  2.11  yields 


E 

z 


1 


JwyoCo 


82A. 


8  z" 


+  k02A. 


(2.13) 


where  ko=  w/poeo  »  the  free  space  propagation  constant. 
The  magnetic  vector  potential,  A(r)  ,  arising  from  a 
volumetric  current  density,  J(r' )  ,  can  be  obtained  from 

/-jk0R 

t(f)T¥r—  dv'  (2-14) 

V' 
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f 


where  R  is  the  distance  from  the  observation  point  (x,y,z) 
to  the  source  point  (x',y',z')  and  is  given  by 


R  =  /(x-xV)z  +  (  y'-'y  1  J  ^  +  (  z  -  z  '  Y  V‘*A3; 

Noting  that  the  free  space  scalar  Green's  function  is 
expressed  as 


G(r.r') 


e  -  j  k  o  R 
4ttR 


(2.16) 


and  that  the  volumetric  current  density  only  has 
a  z  component,  that  is  J(F')szJ  (r1  )  ,  equation  2.13  can 
now  be  written  as 


E 

z 


1  /’[if. 

jwe  o  J  9* 

\/  I 


zT7G(r,r' )  +  k0zG(r,r'  ) 


Jz(r' 


)dV ' 


(2.17) 


This  volume  integral  can  be  further  reduced  to  a  line 
integral  by  making  the  following  assumptions: 

1.  The  current  is  confined  to  the  surface  of  the 
dipole;  this  is  true  for  infinite  conductivity,  and 
approximately  true  for  good  conducting  material; 
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2. 


The  observation  point  is  chosen  to  be  a  point  on 


the  z  axis,  thus  reducing  equation  2.15  to 

R  =  /(T-T'?  +  ‘a 2  (2.18) 

3.  The  current  is  uniform  with  respect  to  <j> 1  , 

a<<A  .  Thus  we  can  write  equation  2.17  in  line  integral 

form  as 


E 


z 


_ 1__ 

j(D£  o 


+ 


L 

2 


6(z,z‘ )  +  k  o  2  G ( z,z‘ ) 


I  (  z ' ) dz 1 


L 

2 


(2.19) 


This  equation  states  that  the  induced  current  can  he  written 


as  a  line  current  located 

a 

distance  a 

from 

the 

observation  point;  note  that 

a 

can  approach 

zero 

but 

cannot  go  entirely  to  zero. 

By 

the 

use  of  the 

Equi valence 

Principle  (Ref  5:104-110), 

Ez 

i  n 

equation  2. 

19  can 

be 

written  as  Es  ,  the  free 

space 

field  radiated  by 

the 

equivalent  current  I(z')  .  Also,  since  a  perfect 
conductor  was  assumed,  there  will  be  no  losses  at  the 
dipole.  Therefore,  the  incident  field,  e1  »  which  created 
the  induced  current  I(z')  .  and  the  scattered 
field,  .resulting  from  I^z')  »  must  be  related  at  the 


19 


dipole  so  that  the  total  electric  field  is  zero,  that  is 


or, 


E1  + 
z 


(2.20) 


(2.21) 


Thus  we  can  write  equation  2.19,  the  Pocklington  Integral 
Equation,  as 


-E 


i 

z 


3 k'/[3P€U-z')+k»2G(z'z') 

L 

2 


I  ( z  1  )  dz  1 


(2.22) 


Equation  2.22  can  be  rewritten  as 


E 


s 

z 


1 

4tt  JU)£  0 


e-jk0R 

Rl 


( 1+ j  k 0 R ) ( 2R2 -3a 2 ) 


(2.23) 


+k„2Rza2 


I  ( z  1  )  d  z 1 
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2.2.2  The  Method  of  Moments.  A  specific  form  of  the 


method  of  moments  is  developed  and  used  herein;  this  form 
being  the  point  matching  technique  using  pulse  expansion 
functions . 

The  dipole  is  divided  into  N  segments,  each  of 

length  ' =  L / N  >  with  each  segment  assumed  to  have  an 
n 

unknown  constant  current.  A  more  accurate  representation 
of  I(z')  can  be  obtained  by  increasing  the  number  of 
segments,  thereby  tending  to  smooth  out  the  stairstep 
approximation;  in  general  N  cannot  be  made  arbitrarily 
large  due  to  finite  word  length  in  the  computer,  leading  to 
numerical  instability.  With  the  above  assumption,  we  can 
write  the  unknown  current  distribution  as 


Hz') 


N 


E 


Fn(z') 


(2.24) 


where  i  is  the  complex  expansion  coefficient  of 

J.  L. 

the  n  segment  and  F  (z')  is  the  pulse  expansion 
function,  defined  as 


Fn(zl)  =< 


1  ;  i f  z '  i s  i n  Az  ' 

n 


0  ;  elsewhere 


(2.25) 
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Also,  for  notational  brevity,  let  us  define 


K  U,z') 


_1 _ 

JU)£  o 


32G(z  ,z'  ) 
3  T2 


ko2G (2,2') 


(2.26) 


Designating  the  center  point  of  the  particular  segment  being 

solved  for  as  zm  ,  the  scattered  field  arising  from 
^  1*1 

the  mc  segment  can  be  written  as 


Ez(zm}  =  f  E  *n  Fn(z'>  K(zm’z'>  dz'  <2*27> 

J  n=l 

L 

‘  2 


If  we  now  use  the  definition  of  F  ( z ' )  ,  we  can  reduce 

n 

equation  2.27  to 


(2.28) 
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and,  if  we  let  f (zml zn> )=  J  K(zm.zn')dzl  we  have 


Az  ' 
n 


Ez<z„>  -  >■  f<VZ>'> 


+  •  •  •  +  IM  f ( z_ , ZM' ) 
N  m  N 


(2.29) 


In  order  to  solve  equation  2.22  at  the  mth  segment,  we 
must  have  the  sum  of  the  scattered  fields  from  all  N 
segments  plus  the  incident  field  at  zm  set  equal  to  zero. 
We  can  therefore  write  N  independent  equations  with  the 
only  unknowns  being  the  N  I  's 
The  resulting  N  independent  equations  being 

-E*(zi)  =  I  if (zj ,Zi ' )  +  I2f ( z ! ,z2 ’ )+  •••  +INf(zi,zN') 

-E^Za)  =  I1f(z2,zi')+  •  •  •  +1  Nf  ( z2 ,  zN '  ) 

(2.30) 

•  •  • 

•  •  • 

-E^(zn)  =  Iif(zN,Zi')+  •  •  •  +INf(zN,zN') 

These  equations  can  be  written  more  conveniently  in  matrix 
notation  as 
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r 

= 

f (Zi  ,Zi  '  )  f (Zi  ,z2 '  )  •  •  •  f (Zi  ,zN'  ) 

II 

-e’(22) 

= 

f (z2  ,Zi  ' )  f (z2  ,Z2  '  ) . 

I2 

• 

• 

• 

-Ez‘V 

•  • 

•  • 

•  • 

f(ZN.Zi')  f(zN,zN') 

_ 

• 

• 

• 

(2.31) 

And  in  abbreviated  form  as 

[V  ]  =  [Z  ][I  ]  (2.32) 

L  L  mnJL  nJ 

where  Vm  =  -Ej(zm)  ’  and  zmn  =  fUm,zn')  the  imPedence 

matrix  which  was  calculated  using  Simpson's  Rule.  As  was 

noted  earlier,  the  purpose  of  this  procedure  is  to  solve  for 

the  i ^  *  s  ;  therefore,  if  we  premultiply  equation  2.32 

by  rz  I" 1  we  obtain 
L  mnJ 

VJ  "  (2'33> 

Since  a  thin  dipole  is  being  used  in  this  analysis,  the 
impedence  matrix  [z  ]  will  be  Toeplitz  (Ref  12:205); 
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I 


therefore,  only  one  row  of  the  matrix  need  be  determined. 

TZ  I”1  can  then  be  determined  using  a  Toeplitz  inversion 
L  mnJ 

routine.  Such  an  inversion  routine,  SUBPROGRAM  TOPIN,  was 
generated  by  modifying  SUBPROGRAM  TPLZ  (Ref  14:579-81),  see 
Appendix  A. 

2.3  Fields  at  the  Antenna 

The  electric  field  present  at  the  receiving  antenna  is 
the  superposition  of  the  field  received  directly  from  the 
leaky  coaxial  cable  and  the  field  that  is  radiated  by  each 
of  the  N  sub-dipoles. 

The  field  received  directly  from  the  leaky  coaxial 
cable  can  be  simply  calculated  using  the  vertical  component 
of  the  electric  field  as  described  by  equation  2.6.  This 
received  field  is  given  by 

E1  =  [E(Bi)  +  E(B2)]  sinUl  (2.34) 

where  <|>  is  the  angle  of  incidence,  see  Figure  2.1. 

The  field  received  from  the  N  sub-dipoles  is  found  by 
examining  each  sub-dipole  separately.  Knowing  that  the 

A 

current  has  only  a  z  component,  and  that  each  segment  is 
small  enough  so  that  it  appears  to  be  a  point  source  of 
amplitude  In(jj-)  ,  solution  of  the  Helmholtz  equation 
yields 
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(2.35) 


A  (r)  -  I  t  e.xJj|-j.k.tI.) 
n  n  N  4irr 


x  L. 

where  r  is  the  distance  from  the  n  n  sub-dipole  to  the 
receiving  antenna.  Now,  by  using  Maxwell's  equations  and 
the  assumption  that  r>>\  ,  we  can  approximate  the  electric 

X  L 

field  at  the  receiving  antenna,  from  the  ntn  segment,  as 


E 


n 


e-jkor 


(2.36) 


where  0  is  the  angle  as  measured  from  the  positive  z-axis, 
and  no  is  the  free  space  characteristic  impedence 
(Ref  5:78-9).  Consequently,  the  electric  field  at  the 
receiving  antenna  from  the  N  sub-dipoles  will  be 


E 


N 


n=  1 


sin(yn) 


(2.37) 


or 
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(2.38) 


N 

E  =  flj-  sin(Q)5I  Ine'jk°rsin(Yn) 

n  =  l 


where  v  is  the  angle  of  incidence  of  the  electric  field 
'  n 

J.  L. 

radiated  from  the  nl  sub-dipole. 

2.4  Vol tage  at  the  Antenna 

The  receiving  antenna  was  arbitrarily  chosen  to  be  0.5 
meters,  or  0.1  wavelengths,  in  length.  Since  this  antenna 
is  located  in  the  far  field  of  the  ambient  and  intruder 
scattered  fields,  there  were  no  non  linear  spatial  effects 
and  the  antenna  was  considered  to  act  as  an  ordinary 
antenna.  The  relative  voltage  was  found  by  comparing  the 
received  field  strengths. 
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III.  RESULTS 


This  chapter  presents  the  theoretically  derived  results 
and  compares  them  with  those  calculated  and  experimentally 
obtained  by  Poirier  (Ref  8). 

3. 1  Theoretically  Derived  Results 

It  was  observed  that  the  electric  field  received 
directly  from  the  coaxial  cable,  equation  2.34,  is 
negligible  when  compared  to  the  electric  field  which  is 
scattered  from  the  intruder-dipole,  equation  2.38.  This  can 
be  explained  by  examining  the  rate  of  decay  of  the  two  field 
components.  Figure  3.1  is  a  plot  of  the  radiated  power 
versus  distance  for  both  the  coaxial  mode  and  the  'Goubau' 
mode;  a  ^  decay  rate  is  also  plotted  as  a  reference. 
Examination  of  this  figure  shows  that  for  p>0.5  meters, 
the  power  begins  to  fall  off  sharply  from  the  rate, 
with  the  rate  of  decay  being  dominated  by  the  exponential 
behavior  of  the  Bessel  function.  The  electric  field  from 
the  coaxial  cable  is  a  non-radiating  field  whose  decay  is 
dominated  by  e"p  ;  where  as,  the  field  scattered  from  the 

dipole  is  a  radiated  field  and  decays  as  —  .  The  field 

p 

received  from  the  coaxial  cable  was  calculated  to  be 
approximately  800  dB  down  from  the  field  scattered  from  the 
dipole,  for  the  case  when  the  receiving  antenna  is  located 
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Figure  3.1  Power  Decay  Rate  for  the  Coaxial  and  'Goubau'  Mode 
The  Straight  Line  is  a  1/p2  Variation. 


12  wavelengths  from  the  coaxial  cable. 

The  variation  in  the  received  power,  as  the  dipole 
moves  within  ±  two  wavelengths  of  the  coaxial  cable  is  shown 
in  Figure  3.2.  This  figure  shows  the  relative  power 
received  at  the  antenna  for  various  intruder  distances  from 
the  cable;  where  the  received  power  is  normalized  to  the 
maximum  power  level,  the  power  received  when  the  intruder  is 
at  the  cable.  This  curve  shows  that  the  received  power 
falls  off  drastically  with  distance  from  the  cable;  and, 
that  intruder  detection  would  be  possible  with  a  very  narrow 
detection  region.  Figure  3.3  shows  the  received  voltage, 
again  normalized  to  the  value  received  when  the  intruder  is 
at  the  cable,  for  the  probable  zone  of  detection,  ±  10  feet 
or  +0.6  wavelengths . 

As  noted  above,  the  theoretical  calculations  indicate 
the  electric  field  received  from  the  coaxial  cable  to  be 
negligible;  however,  Poirier  reports  that  the  ambient  field, 
or  the  field  received  from  the  coaxial  cable  alone,  is  of 
the  same  order  magnitude  as  the  maximum  field  scattered  from 
the  intruder-dipole  (Ref  9).  For  this  reason,  a  separate 
plot,  Figure  3.4,  was  made  wherein  the  field  received  from 
the  coaxial  cable  was  arbitrarily  chosen  to  be  37%  of  the 
field  scattered  from  the  dipole  when  it  is  located  at  the 
coaxial  cable,  so  that  a  comparison  to  Poirier's  results 
could  be  made.  In  Figure  3.4,  the  dip  occurs  when  the 
ambient  field  and  the  intruder  scattered  field  are  out  of 
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Figure  3.3  Normalized  Received  Voltage  Variation  as  an 
Intruder  Moves  Along  a  Radial  Path.  (Note:  The  "glitch" 
near  the  cable  is  not  really  there  and  is  due  to  problems  in 
the  plot  routine  when  the  dipole  is  directly  over  the 
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D I  STANCE  FROM  CRBLE 


phase  and  thus  subtract;  the  peak  occurs  when  they  are  in 
phase  and  the  fields  add. 

3 . 2  Review  of  Poirier  Report  (Ref  8) 

Poirier  has  presented  a  method  for  explaining  the 
received  field  variations  due  to  the  interaction  between  an 
intruder  and  the  fields  produced  by  the  leaky  coaxial  cable, 
experimental  measurements  are  also  incorporated. 

Poirier's  approach  makes  use  of  the  following 
assumptions/conditions: 

1.  The  coaxial  cable  is  located  on  the  surface  of  the 
earth  ; 

2.  The  earth-air  interface  does  not  significantly 
modify  the  relative  system  response  (another  way  of  saying 
free  space  analysi s ) ; 

3.  The  electric  field  resulting  from  the  leaky  coaxial 
cable  is  a  combination  of  the  coaxial  mode  and  the  'Goubau' 
mode ; 

4.  The  intruder  is  modeled  as  a  point  source  located 
0.4  meters  above  the  ground;  and, 

5.  The  frequency  of  operation  is  75  MHz. 

Poirier  then  goes  on  to  conclude  that  the  power  density 

falls  off  as  ^-z-  for  all  n's  of  interest;  this  conclusion 

P 

is  based  on  a  comparison  of  the  coaxial  and  'Goubau'  mode 
power  densities  with  a  ^  decay  rate  (similar  to  Figure 
3.1),  and  data  reported  by  Cree  and  Giles  (Ref  3)  indicating 
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that  the  power  decay  rate  is  ~r  for  all  p 

Finally,  Poirier  makes  use  of  measurements  indicating 
that  the  ambient  electric  field,  the  field  received  when  the 
intruder  is  absent,  is  of  the  same  order  magnitude  as  the 
field  that  is  scattered  from  the  intruder  (Ref  9).  So,  by 
letting  the  total  received  electric  field  be  the  sum  of  the 
ambient  field  and  the  intruder  scattered  field,  and  by 
normalizing  to  the  ambient  field,  he  obtained  the  power 
variation  shown  in  Figure  3.5,  which  agrees  closely  with 
experimental  measurements. 

3.3  Comparison  of  Results 

Based  on  the  results  of  sections  3.1  and  3.2,  the 
inescapable  conclusion  is  that  a  "free  space"  analysis  is 
incompatible  with  experimental  results  which  dictate  a  power 
decay  rate  that  is  primarily  determined  by  .  This  is 

further  pointed  out  by  realizing  that  for  the  point  source 
intruder  model,  at  the  point  of  closest  approach,  0.4 
meters.  Figure  3.1  shows  that  this  is  the  distance  where  the 
power  decay  rate  begins  to  no  longer  follow  ~z-  .  One  can 

only  conclude,  that  the  presence  of  the  earth-air  interface 
somehow  enables  a  form  of  a  bound  wave  to  propagate; 
possibly  accounting  for  the  ^  variation  instead  of 
the  Ki(p)  variation  for  the  radial  electric  field 
component. 

Comparison  of  Figures  3.4  and  3.5  again  show  that  the 
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! 

chosen  form  of  the  field  structure  is  not  correct  since  the 
received  power,  of  Figure  3.4,  is  just  the  ambient  power, 
except  when  the  intruder  is  on  top  of  the  coaxial  cable. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


4. 1  Cone! usi ons 

The  intent  of  this  investigation  was  to  determine  the 
free  space  interaction  between  a  dipole  and  the  EM  fields 
created  by  a  leaky  coaxial  cable.  In  order  to  accomplish 
this,  an  algorithm  for  determining  the  effects  of  a  non¬ 
plane  wave  incident  upon  a  dipole  was  developed.  It  was 
then  found  that  by  exciting  a  current  on  a  dipole,  the  far 
field  radiation  would  be  greatly  enhanced  by  converting  the 
evanescent  field  of  the  coaxial  cable  into  a  radiated  field 
from  the  dipole.  However,  since  this  study  was  undertaken 
in  order  to  analyze  the  interaction  of  an  intruder  with  the 
EM  fields  created  by  a  leaky  coaxial  cable,  the  following 
conclusions  are  made: 

1.  Intruder  detection  is  possible.  The  zone  of 
detection  is  just  smaller  in  the  free  space  case  (which  may 
be  better). 

2.  The  effects  of  the  earth-air  interface  cannot  be 
ignored  and  still  obtain  exact  results.  The  general 
character  of  the  response  is  similar  but  the  ambient  levels 
very  significantly. 

3.  The  form  of  the  EM  field  is  apparently  affected  by 
the  earth-air  interface. 
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4.2  Recommendations 


Based  on  the  assumptions  stated  earlier,  and  the 
results  presented  herein,  the  following  recommendations  are 
made.  This  problem  should  be  pursued  using  the  approach 
presented  by  Van  der  Pol  (Ref  15),  with  a  leaky  coaxial 
cable  as  the  source,  in  order  to  determine  the  proper  form 
of  the  external  EM  fields.  Another  possible  reference, 
using  a  leaky  coaxial  cable  as  the  source  (Ref  7),  could  be 
used  in  conjunction  with  Van  der  Pol  to  solve  the  problem  of 
the  cable  on  the  earth  instead  of  buried,  as  done  in 
Reference  7.  Once  the  fields  have  been  properly  determined, 
the  computer  program  developed  as  part  of  this  study. 
Appendix  A,  can  be  used  to  determine  the  interaction  with  an 
intruder.  Once  the  proper  form  of  the  EM  fields  are  known, 
it  should  then  be  clear  what  effects  the  changing 
environmental  conditions  would  have;  this  is  an  extremely 
important  aspect  of  system  performance  which  needs  further 
i nvesti gati on . 

After  the  above  has  been  accomplished,  one  could  then 
begin  to  determine  a  more  realistic  model  for  the  intruder. 
The  modeling  of  biological  bodies  for  EM  radiation  has  been 
undertaken  by  Chen,  Livesay  and  Guru  (Ref  2),  and  by  Livesay 
and  Chen  (Ref  6) . 


39 


Bib! iography 


Abramowitz,  Milton  and  Irene  A.  Stegun.  (1964)  Handbook 
of  Mathematical  Functions.  U.  S.  Government  Printing 
Office,  Washington,  D.C.  (National  Bureau  of  Standards 
Applied  Mathematics,  Series  55). 

Chen,  Kun-Mu,  Donald  E.  Livesay,  and  B.S.  Guru. 
"Induced  Current  in  and  Scattered  Field  from  a  Finite 
Cylinder  with  Arbitrary  Conductivity  and  Permittivity," 
IEEE  Transactions  on  Antennas  and  Propagation.  Vol  AP- 


Prac 

ectr 


neti  c 


RADC-TR-79-17 ,  in 
house  report.  Hanscom  AFB,  Massachusetts:  Deputy  for 
Electronic  Technology,  Hanscom  AFB,  January  1979.  (AD 
A072  818) 


40 


11.  - et  al .  VHF  Intrusion  Detection:  A  Technique 

for  Parked  Aircraft.  RADC-TR-77-384 ,  in  house  report. 
Hanscom  AFB,  Massachusetts:  Deputy  for  Electronic 

Technology,  Hanscom  AFB,  November  1977.  (AD  A051  144) 

12.  Preis,  Douglas  H.  "The  Toeplitz  Matrix:  Its  Occurrence 
in  Antenna  Problems  and  a  Rapid  Inversion  Algorithm," 
IEEE  Transactions  on  Antennas  and  Propagation,  Vol  AP- 
20:  204-6  (March  1972). 

13.  Rao,  K.V.N.  e t  a  I .  "Resource  Protection  by  Segmented 
Leaky  Coax  Cables,"  Proceedings  1979  Carnahan  Conference 
on  Crime  Countermeasures:  133-7.  Lexington,  Kentucky: 
University  of  Kentucky,  May  16-18,  1979. 

14.  Stutzman,  Warren  L.  and  Gary  Thiele.  Antenna  Theory  and 
Des i gn .  New  York:  John  Wiley  and  Sons,  1981. 

15.  Van  der  Pol,  B.  "On  Discontinuous  Electromagnetic  Waves 

and  the  Occurrence  of  a  Surface  Wave,"  IRE  Transactions 
on  Antennas  and  Propagation,  Vol  AP-41  2’88-93  (July 

1956). 

16.  Von  Hippel,  Arthur  R.  Dielectric  Materials  and 
Appl i cations .  New  York:  John  Wiley  and  Sons,  1958. 


onrtoooftnonoonnoonooooonnonnnnononnonnnooftonoooonooo 


Appendix  A 


Developed  Computer  Program 


PROGRAM  SkIRPS 

***»*$«»*$»*e*$*»»**»*********** 

*  THIS  PROGRAM  CU-PUTci  T  Hr:  FLEC  TKUMAGNc  T  I C  FIELDS 

*  ARISING  FROM  THt-  INPmACTUN  fIF  A  LEAKY  COAXIAL  CABLE 

*  AM)  THUSF  S  C  A  T  I :  <  f  i<  FROM  A  DIPOLE. 

*  THE  FIELDS  SCATTtRf  0  FR-J-  THE  OIPHLE  ARE  CALCULATE!) 

*  USING  T Hl  METHOD  OF  '•U',tNrS  kITH  PULSE  FUNCTIONS. 

*  NOTE:  ALL  01  STANCES  AKF  NORMALIZED  TO  UNITS  OF 

*  WAVELENGTHS,  AND  THUS  BECOME  0 1 ME NS  I ONL E SS  NUMBERS. 

*  LEGEND 

*  AA  »  LOwER  L I  * l T  OF  K-TH  ELEMfcNT 

*  AB<  )  «  MIDDLE  OF  THE  K-TH  ELEMENT 

*  ABAH  »  DIPOLE  RADIUS 

*  ADIST  «  DISTANCE  FROM  COAX  TO  CENTER  OF  K-TH  ELEMENT 

*  BH  *  UPPER  LIMIT  OF  K-TH  ELEMENT 

*  3DIST  -  DISTANCE  FROM  RECEIVE  ANTENNA  TO 

*  CENTER  OF  K-TH  ELEMENT 

*  BETA  -  PROPAGATION  CONSTANT 

*  bETAl  «  COAXIAL  MOOE  PROPAGATION  CONSTANT 

*  «  1.72  RAOS/M 

*  SETA2  «  •GOUBAU*  MODE  PROPAGATION  CONSTANT 

*  »  1.32  RADS/* 

*  COIST  «  OISTANCt  FROM  RECEIVE  ANTENNA  TO  COAX 

*  C UR (  I  «  INDUCED  CURRENTS 

*  Cl  »  RATIO  OF  E  FIELO  FOR  COAX IAC/GOUBAU  MOOES 

*  ODISTC  I  -  DISTANCE  From  COAX  TO  INTRUDER 

*  ED  1ST  -  DISTANCE  FROM  COAX  TO  CENTER  OF  KA-TH 

*  ELEMENT  OF  RECEIVF  ANTENNA 

*  ELBAR  -  DIPOLE  LENGTH 

*  E  RE  C(  I  »  6  FIELO  AT  THE  RECEIVE  ANTENNA 

*  ESCAT1  »  -  E  FIELD  SCATTERED  FROM  THE  K-TH  ELEMENT 

*  K  *  ELEMENT  NUMBER  FOR  SOURCE  ELEMENT 

*  N  •  NUMBER  OF  EXPANSION  FUNCTIONS 

*  PHASE  »  PHASE  DIFFERENCE  BETwEEN  COAXIAL  £  GOUBAU  MOOES 

*  PGW(  I  «  RECEIVED  PU*fcR 

*  R8AR  »  OISTANCE  FRU“  SOURCE  TO  OBSERVATION  POINT 

*  TAUI  I  -  FIRST  ROW  OR  COLUMN  OF  THE  I MPfc  DENCE  MATRIX 

*  U  •  DISTANCE  VARIABLE 

*  VIN(  »  •  INCIDENT  E  FIELO,  VERTICAL  COMPONENT 

*  VCLTSI  »  -  ANTENNA  VOLTAGE 

*  zinvi  i  ■  inverted  impedence  matrix 

*  20  -  CHARACTERISTIC  IMPEDENCE  (FREE  SPACE) 

*  20BS  -  OBSERVATION  POINT 


COMMON  A8AR ,  BETA,  Z03S,  ZJ 

COMPLEX  A  1 1 00  )  ,  All  100)  ,  CURUOO),  ERECIIO),  ESCAT(IOO),  TAU(IOO), 
1  V  INI  100 )  ,  VCLTSU2),  Z  INVI  100,  100) 

COMPLEX  ADD,  AI.MONE,  CONE,  CZERO*  FVAL,  VMAX,  2IJ,  Z J 
INTEGER  COLMN,  ROW 
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REAL  A3I100I.  00151(42).  P0KI42) 

REAL  <1 

DATA  AIMuNt/IO. 0,1.01/,  C  0  N  6  /  (  L.o.O.OI/t  Cl  E  RO  /  (  0. 0  »  0. 0 »  / 

ABAR-U.i.08 

CL*  1.0 

COIST-12.0 

ELSAR-G.  3t>8 

idone-o 

INVERT-O 
N*  100 
PHASE-0.0 

P  I  *  3.  14  1  5  92  68  38  89  3 
VMAX-C2  fc  '<  0 
20-376.7 
2J-AIM0NE 
BETA-2. 0*PI 

8ETA1  »  1.72*8.0 
BET A2  -  1.32*8.0 
0EL-EL6AR/N 
20BS-ELBAR/ 12 .0*N) 

00  90  L--21..20.1 
L V  -  L*  22 

DOISTILVt  -  0.08  ♦  L/ 10. 0 
00  20  K  ■  1  .  N 

IF  I  1 OONE •  NF • 0 1  GU  TO  18 

DETERMINE  LOWER  LIMIT  OF  K-Th  ELEMENT 

AA-ELBAR*(K-1.0I/N 

OETERMINt  UPPER  LIMIT  OF  K-TH  ELEMENT 
flB-EL  8AR*K /N 

OETERMINE  CENTER  OF  K-TH  ELEMENT 
A8IKJ-IAA.8B1/2.0 

COMPUTATION  OF  IMPEDENCE  MATRIX 

PERFORM  INTEGRATION  USING  SIMPSONS  RULE.  REF  "CALCULUS  AND 

ANALYTIC  GfcOMETR Y.2E0".  BY  JOHN  F  RANOOLF,  PG  427.  EQ  5. 

NST  EPS-80 

NN-NST  EP  S/2 

OINC-IBB-AAI/NSTEPS 

AOO-C2ER  0 

U-AA 

DO  10  J-l.NN 
CALL  FUNCTIU.FVALI 
AOO«AOO*2.0*FVAL 
U-U+OINC 

CALL  FUNCTIU.FVALI 
A00»A0D.4.0*FVAL 
U-U*0 I NC 
CONTINUE 

AOO  IN  LAST  TERM  ANO  CORRECT  FIRST  TERM 

U-AA 

CALL  FUNCTIU.FVALI 
AOO-AOO— F VAL 
U-BB 

CALL  FUNCTIU.FVALI 
AOO-AOO*FVAL 
21 J- AO 0*0  INC 7 3.0 


T*um«n  j*zo/<  z  j*b.o*p  i*pi  » 

15  CONTINUE 

c 

C  COMPUTATION  CP  I NC  JOt  NT  VERTICAL  t  FIELO 

c 

HQ  =  Si” T  C  BP  T  A 1  ** 2 . 0— S(  rA$#2.0> 

AD1ST  »  SURTCOoliT (LV» v*2.0»  Art(K)**2.0> 

A.<G*MU*AQI'iT 

CALL  PflUdE S ( AwG ,K  l ) 

SINE  A*  AIM*  I  Z  A  0  l  b  T 

VIMK  )»(-2.3*Bl  TA|*K  L/C  HG*P  l )  )*SI  NE  A 
HQ“S0R  T(8fcTA<?®*?.0-'3(  T A**2.0> 

ARG=HO*ADI S  T 

CALL  «ODbESCAPG,Kll 

VlNCK|«VlN<K)-2.*3ETA2*Kl*SINEA*CEXPC-ZJ*PHASE)/ChO«Pt*Cl) 

20  CONTINUE 
I  DONE “l 
C 

C  COMPUTATION  OF  2— INVERSE 

c 

IF  C INVERT. NE.O)  GO  TQ  30 

INV£RT*1 

NZ»N 

CALL  TOPINVITAUtNZiA.Al.ZINV) 

30  CONTINUE 

C 

C  COMPUTATION  OF  INDUCED  CURRENT 

c 

DO  50  RO»»l,N 

cur c  row ) *c  zero 

00  50  CQLMN*= 1 ,N 

CUR  CRON »-CUR( RO  W) »Z l  N V C RQH ,  COL KN ) *  V I N CC QL HN  ) 

50  CONTINUE 

50  CONTINUE 

c 

C  DETERMINE  THE  VERTICAL  E  FIELD  AT  ANTENNA  DUE  TO  COAX 

C 

00  60  XA-1,10 

ANT«K A/10.0-0.05 
PHl*ATANCANT/COIST) 

EQlST«SORTCANT**2.0»C0IST**2.O) 

H0*5Q9TCeE  TAl**2.0-dETA**2.0> 

ARC  l  *HO»E 0  1ST 

CALL  MaOBESCARGl.Kl) 

ERECCXAI-»-2.0*3ETAl*Kl/ CHO*P I  >  >*SINCPHI ) 

H0» SOR T ( BE TA2**2.0-3ET A* *2.0) 

ARG1«H0*EQIST 

CALL  MOOBESCAHGl.Kl) 

ERtCCKAl  «  C-2.G*HETA2*Kl/CH0*PI*Cl) >*SINCPHI > *CEXP C -Z J*PHASE ) 
l  ♦  ERECIKA) 

60  CONTINUE 

C 

C  COMPUTATION  OF  THE  SCATTEREO  VERTICAL  E  FIELD  AT  THE  ANTENNA 

C  FRO**  THE  K-Th  ELEMENT  AND  THE  RESULTING  RECEIVED  E  FIELO 

c 

VOLTSCLVI  -  CZERO 
00  to  K-l.N 
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00  70  KA-l.LO 

A  N’T  *KA/IU.0-0.G5 
01  f  =  AHS  (  Att  CK  l-A'ill 

B0IST<S)><T(')lF«*7.l>*IC0IST«iJ01ST(LV))**2.0l 
GAPBA  =  AS1MA!<  (  <  )  /-*  0  l  S  T  ) 

TriETA  =  GAB“A,Pl/i'.0 

FSCAT|K>=2  J«CJ-(»K>«,'FL*/U»SIN<THFTA»*CfXP(-Z  J*BETA*BOISH/ 
i 

CHI » AT  AN t  D I f / I  CO  I S  T ♦ UQ  I  ST l L  V  1 1  I 
ERECC  <A  J  =Ek£C (Ki)-L  SC AT(K  )»  CUS ( CH I  ) 

VULTS(LV)  »  YOLTS(LV)  ♦  0 .  0  l*h.R  t  C  t  K  A ) 

70  CONTINUE 

80  CONTINUE 

IF  (CASS  I  VOLT  S  C  LV  J I  .GT.  CABS(VPAX1I  VPAX  *  VOLTS(LV) 

90  CONTINUE 

DO  100  LP “ 1 , A  2 

POp(LP)  -  l0.U*L0GLC{(CArtSIV0LTS(LP>l/CABSlVPAXll**2.01 
100  CONTINUE 

PRINT* (T5,A,T20,AI*. ’INTRUDER**’ RECEIVED* 

PRINT* (T5  *A,T21,A) *  , 'DISTANCE*  ,*PJ..ER • 

PRINT*  (T7«F5.2,T20,F7.2»  *»  (001  ST  (NO)  ,PO»I P0>  , PD- 1,42) 

STOP 
f  NO 

SUBROUTINE  TOPI NV( T AU.N2 , A , A l . 2 1 N V  » 

* 

*  THIS  SUBROUTINE  DETERMINES  THE  INVERSE  OF  A  T0EPLIT2 

*  MATRIX  WHEN  GIVEN  EITHER  THE  FIRST  ROh  OR  COLURN  OF  THE 

*  T0EPUT2  PATRIX. 

* 

*  NOTE:  THIS  ROUTINE  IS  BORE  EFFICIENT  than  A  GENERAL 

*  MATRIX  INVERSION  ROUTINE  BECAUSE  IT  PARES  USE  OF  THE 

*  SPECIAL  PROPERTIES  OF  A  T0EPLIT2  MATRIX. 

* 

*  THIS  SUBROUTINE  IS  A  MODIFICATION  OF  SUBROUTINE  TPL2 

«  FROP  "ANTENNA  THEORY  AND  DESIGN" 

*  8 Y  STUT/PAN  ANO  THIELE.  PAGE  579  -  81. 

* 

*  LEGENO 

* 

♦  A » A 1  -  SCRATCH  AREA  VECTURS  OF  LENGTH  N2 

•  IER  *  ERROR  CUOt*  FUUAL  TU  2EkO  IE  PATRIX  IS  INVERTIBLE 

♦  N2  OROER  OF  THE  T0EPLIT2  "ATRIX 

*  TAU  *  FIRST  KOb  UR  COLUPN  OF  THE  TUEPLIT2  PATRIX 

*  2  I N V  «  THE  INVERTED  TOEPL I T  2  PATRIX 

* 

**«»»»**#«•«**«****<»«««********* 

CQPNON  ABAR,  BETA,  208S,  2J 

COPPLEX  A  (  N  2  I  *  AKN21,  TAUIN2I.  2INVIN2.N2) 

COPPLEX  ALPOA,  ALPHA,  Cl,  C2,  COEF,  CONE,  C2ER0,  FAC,  TAUl 
OATA  CONE/I  1.0,0.012,  C2ER0/ 1 0, 0, 0. 01 / 

N-N2-1 
IER-0 

TAUl'TAUIl) 

00  2000  1 1 •!  ,N 


2000  TAU(III*TAU(I l«l)/TAUl 

ALN0A»C0NE-TAU t  1  )*T AUl l) 

AC1)*-TAU<1) 

1*2 

1  KK*l-L 

ALPHA*CZERO 
00  2  M=1,KK 
LL*I-(1 

2  ALPHA*ALPHA*A(M)*TAU(LL> 

ALPHA*-! ALPHA*TAU( I ) ) 

IF  (CABS! ALPHA) .tO. 0.00)  GO  TO  IS 

C0EF=ALPHA/ ALMDA 

ALSO A* ALMOA-COEF* ALPHA 

00  3  J=1 »KK 

L*I-J 

3  Al(J)-A( J)»COEF*A(L) 

DO  7  J  *  l »  KK 

7  A(J)*AMJI 

A ( I l-COEF 

IF  (I.GE.N)  CO  TO  5 

4  1-1*1 
CO  TO  1 

5  NH*(N7»l)/2 
FAC*ALHOA*TAUL 
NP  *N2 ♦ l 

00  SI  1=1 ,NH 

IF  (  I.NE.l )  GO  TO  52 

A1(1I*C0SE/FAC 

00  53  J=2,NZ 

All J)=A( J-ll/FAC 

53  CONTINUE 
GO  TO  54 

52  CONTINUF 

JH= l -I 
C1*A( JH) 

NNP I =NP— I 
C2* A(NNPI ) 

00  55  JJ«l,N 
J*NP-J J 
INP J=NP-J 
JL  * J— 1 

AIU)=AI  (  JL  I  *(CI«A(  JL)-C2*A(  INP  J)  l/FAC 
55  CONTINUE 

Aim=A(  i-n/FAC 

54  CONTINUE 

00  57  J*l,NZ 
2INV( I , J l*AL( J ) 

KNPJ*NP-J 

KNPI-NP-I 

2 1 NV ( AN  P I «  J ) • A i ( KNP  J ) 

57  CONTINUE 

51  CONTINUE 

RETURN 

15  P  R I  NT • E  RROR  HAS  OCCURREQ.  HATRIX  IS  NOT  STRONGLY  NUNS  INC* 
I E  R  *  I 
RETURN 
ENO 
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SUBROUTINE  FUNCTIU.FVALI 

t*»t»»#<.0(,45»4fs«e**e*********** 

* 

*  THIS  SUBROUTINE  CUMPUTtS  THE  VALUE.  OF  THE  FUNCTION  BE  INC 

*  INTEGRATED  IN  EQUATION  2.23. 

* 

«««*««*  *  <■  *  «  * 

COM*uN  AHAR,  BETA,  EOaS,  2J 
COMPLEX  Al,  FVAL,  2J 
U  TU2=U-/U8S 

R8AR*SJrTIA8AR**2.0*UT(I/<'#2.0I 

Al  =  (2.GRKBAK**2.C-3.0*AflAR**2.0J*(  1.0«/J*8ETA*RaARI 
A2»< B£TA*R8AR *Ad  AR »**2.0 

FVAL-<CEXP«-/J*3£T A*RHAR I*IA1*A2 )>/(kBAR**5.0) 

RETURN 

END 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


c 


1000 

C 


SUBROUTINE  MOOBE  S (  X ,K  1 ) 

********************************* 

♦  * 

♦  THIS  SUBROUTINE  DETERMINES  THE  VALUE  OF  THE  MODIFIED  • 

♦  dESSEL  FUNCTION  OF  THE  SECONJ  KIND,  ORDER  ONE,  GIVEN  • 

*  A  REAL  ARGUMENT  * 

•  * 

*  legend  * 

*  * 

♦  II  *  MQOIFIEO  BESSEL  FCTN  OF  THE  FIRST  KINO,  ORDER  ONE  • 

♦  K1  MODIFIED  BESSEL  FCTN  OF  THE  SECOND  KINO,  OROER  ONE  • 

♦  X  REAL  ARGUMENT  * 

»  ^  4  * 
*♦*********♦*♦******•****♦*••*•♦* 

REAL  11, Kl 

IF  (X.GE. 2.0)  CO  TO  1000 

REFERENCES  ABRAM 0 RITE,  60S  9.8.3  E  9.8.7,  PCS  378-9. 

T'X/3.75 

II «X *1  0  •  5  +  0. 8 789059 4*T**2. 0*0.8 1 498  8f>9*T** 8.0*0. 1508  49  34*T **6 . 0* 

l0.02N5d733RF**B.O*0.0030l532*T**lQ.ON0.00032«ill*T**12.0) 

Kl»(l.0/XJ*(X*L0G(X/2.G)*I !♦ 1.0*0. IB A A  31 A A* IX/ 2.0) ♦♦2,0-0. 6 72 78579 
l«(X/2. 01**4. 0-0. IH15689/*(X/2.0I**6.0-0.01919402*IX/2.0>**8.0- 
20. 00110404* 1X72.01**10.0-0.00004686*  (X/2.0  1**12.01 
RE  TURN 
CUNT INUE 

REFERENCES  A8RAMOWIT2,  EOS  9.8.8,  PCS  379. 

Y»EXPIX)*SORT(X» 

YY* 1.0/ Y 

Kl»<  YY  J*( 1 .253 31 41 4*0. 2349  8619* ( 2. 0/X I— 0.0  3655620*  12 ,0/X  1**2 .0* 
10. 0 150 A 26 8* (2. 0/X 1**3 .0-0 . 00 780 35 i* I2.0/X I* *4. 0*0 .003 25614* I  2 .0/X ) 
2**5 . 0—0 . 00068245*1 2.0/XI**6.0I 
RETURN 
ENO 
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